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cn Abstract 

Force extension curves (FECs), which quantify the response of a variety of biomolecules subject to 
mechanical force (/), are often quantitatively fit using worm- like chain (WLC) or freely-jointed chain 
(FJC) models. These models predict that the chain extension, x, normalized by the contour length 
increases linearly at small / and at high forces scale as x ~ (1 — / a ) where a= 0.5 for WLC and unity 
for FJC. In contrast, experiments on ssDNA show that over a range of / and ionic concentration, x 

o 

scales as x ~ In /, which cannot be explained using WLC or FJC models. Using theory and simulations 

^ we show that this unusual behavior in FEC in ssDNA is due to sequence-independent polyelectrolyte 
> 

effects. We show that the x ~ In/ arises because in the absence of force the tangent correlation 

m 

function, quantifying chain persistence, decays algebraically on length scales on the order of the Debye 
length. Our theory, which is most appropriate for monovalent salts, quantitatively fits the experimental 

o 

CN) data and further predicts that such a regime is not discernible in double stranded DNA. 
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Introduction: 

The response of double stranded DNA, RNA, proteins and polysaccharides to mechanical 
force (/) has provided a microscopic basis for describing their elasticity^'. A number of studies 
have established that the measured force-extension curves (FECs) can be nearly quantitatively 
reproduced using standard worm-like chain (WLC) modePor freely-jointed chain (FJC) modeP. 
However, recent single molecule experiments on single stranded DNA (ssDNA) showed that the 
measured FECs exhibit unexpected behavior that cannot be described using WLC or FJC model. 
Surprisingly, over a wide range of stretching forces / and at low ionic strengths J, the chain 
extension normalized by the contour length of the chain, L, scales as x ~ In yMUl 

The unusual behavior found for a range of / values in ssDNA implies that the standard 
polymer models that take only the elasticity of the chain are not adequate. Here, we provide 
a theoretical explanation of the x ~ In / behavior using analytical calculations and simulations 
by treating ssDNA as a polyelectrolyte chain. We show that the tangent correlation function 
from which the chain persistence length is extracted decays as a power law on scales a < s < 
(the Debye length). Such a behavior, which was not previously noticed, is a consequence of 
polyelectrolyte (PE) effects in ssDNA and gives rise to the x ~ In / behavior. We confirm the 
theoretical predictions using simulations, which also show excellent agreement with experimental 
data for ssDNA^ for I from ImM to 3M. Our theory, not only produces quantitative agreement 
with experiments, but also establishes the importance of PE effects in predicting the response 
of ssDNA to /. 

Theory for stretching ssDNA including polyelectrolyte effects. 

Measured FECs for a broad class of biopolymers can be understood in terms of an interplay 
of the tensile length £ t = fcgT// (fcgT is the thermal energy), the persistence length £ p and the 
monomer length a. The resulting theory allows us to determine a priori the precise polymer 
model that best describes the measured FECs for a specific macromolecule^. For a discrete 
semiflexible chain, the tangent correlation at / = decays as 

(cos0( S )) = exp(-£) (1) 

where 8(s) is the angle between two tangent vectors that are separated by length s along the 
chain (Fig. KLa). At very high forces, / > Th = ckBT^ p /a 2 , {cos 9) is independent of s and can 
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be approximated as^ 

6 



(cos#)(/) ~ exp y—c 

where c is a model-dependent constant on the order of 4 for non-charged polymers. Thus, the 
macromolecule behaves as a FJC with x ~ 1 — £j/a ~ 1 — When ksT/C, p < f < J-jj, due 
to the interplay of £ p and £ t , the macromolecule can be modeled by a sequence of independent 
segments whose length A A is force-dependent. In this force regime, the tangent correlation is 
given by^ 

(cose>(/)wexpf-c^-V (3) 



Thus, A_r- = yjc^t^p and x « 1 — ^/A A - ~ 1 — J -1 / 2 , which is the well-known WLC behavioral 
Exponential decay of (cos(6(s)) in uncharged polymers (Eq. 0) results from the sequential 
transmission of interactions along the chain backbone . However, in polyelectrolytes such as 
ssDNA, direct repulsive interactions between any pair of monomers within a scale ~ Ad could 
alter Eq. ([!]). It has already been shown that deviations from Eq. Q could occur even in stiff 
chains without intra molecular interactions^. The PE effects can profoundly change the nature 
of FECs. Because the anticipated decay of {cos 6) in PEs, which is slower than in uncharged 
polymers affects FEC, we propose that the intrinsic (/ = 0) tangent correlation on scales s < Ad 
in polyelectrolytes should decay as 

(cos0(s)> ^Cf-Wor a<s< X D , (4) 



, a / 

were C and 7 are functions of a/ Ad oc I 1 / 2 . Note that, in the complete screening limit Xd/cl <C 1, 
the correlation equation must reduce to the exponential form, thus (cos#(s = a)) = C = 
exp (— a/£p), with £° being the bare persistence length; C is, thus, the lower limit for C. 

Just as for macromolecules with short-ranged interactions between the monomers, the FEC 
for PEs should quantitatively match the prediction of the FJC model as long as / > JF^ . We 
determine by equating Eqs. ^ and (jij at s — a yielding 

t pe _ k BT 

In the complete screening limit J-[ E in Eq. ^ reduces to Th for uncharged polymers described 
above. 
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In the regime where the force is high enough such that every segment on the order of Ad 
is independent of each other, the chain can be modeled as a FJC with force-dependent Xk 
segments, with the tangent correlation due to the tensile energy being given in Eq. (|3j). By 
equating the rhs's of Eqs. Q and ^ at s = Xk we have the solution for A#-(/), 

Xk = 7 ™ I t\ f ' ^ 

where C* = C~ 1 ^' and W(z), the Lambert Omega function^, is the inverse of z = We w . The 
equation for the FEC for polyelectrolytes becomes 

,«l-f=l-W*™£), (7) 
X K c V 7 fa J 

or equivalently, 



/ « C*e" c/7 — - P ^ X 



i — x 

Thus, in the regime where x is not too close to 1, the exp factor dominates in Eq. (J8|, 

which implies that / ~ exp(x), or equivalently x ~ In /. The condition for observing the unusual 
behavior is c/7 ^> 1. We show using simulations (see below) that when I/7 ~ (Ad/a) ' 36 , this 
condition is satisfied for Ad ^> a or at low ionic strengths, just as observed in experiments^^!. 
Simulations: 

In order to validate the theoretical predictions we first performed Monte Carlo (MC) simulations 
of a PE chain with bending rigidity and screened electrostatic repulsions between the monomers, 
model that is appropriate for ssDNA. The chain Hamiltonian is given by 

1=0 8=0 J 1=0 

j=i+l 

where N is the number of monomers, bj is the z-th bond vector (Fig. [IJi), is the 
distance between beads i and j, A# = 0.7 nm is the Bjerrum length^, and q = 
a exp (0.0338 + 1.36/ 2/5 ) nm 1 with / measured in molarity units (M), and q is the effective 
charge number per monomer in monovalent solutions^. Here a = 0.55 nm and £° = 0.65 nm, 
which are similar to the parameters for ssDNA/RNA'^^HSEQ] D eD y e length is related to / 
(in M) as Ad = 1 / v / 8 71 ~Ad0.602J in monovalent solutions^. 



Sampling of chain conformations is done using the Metropolis scheme with the crankshaft 
and pivot moves^M^. Starting from an initial conformation, the chain is equilibrated for at 
least 10 6 MC steps before the actual sampling starts. The long range nature of the electrostatic 
interactions, even though dampened by the Debye-Huckel factor, is computationally demanding. 
To circumvent this problem, we use a matrix to store all the pair-wise interactions and only 
update the pairs whose distances are changed by the MC moves. To expedite the computations 
we introduce a cut-off distance of 20 Ad, which is large enough that it does not introduce any 
errors. We checked the accuracy by comparing FECs and the tangent correlation data at / = 1 
and 20mM with cut-offs of 10 and 20Ad and no cut-off, and found virtually no discrepancies. 
With these two enhancements the sampling can be sped up by a factor of at least 5 for large N. 

We calculated the tangent correlation for the model polyelectrolyte for I ranging from 1 to 
500mM. To reduce finite size effects, we simulated for each / several replicas of the chain with 
different lengths N, starting from 100 and doubling it every time, until the results for (cos(8(s)) 
are independent of chain length (N > 800) (see Fig. [TJd for data at / = ImM). The smaller I is 
to (the longer Xd) the stronger are the finite size effects. Hence, at low ionic strengths N has to 
be sufficiently large to obtain converged results. From Fig. [Tj^, it is clear that (cos#(s)) decays 
exponentially at large s as commonly accepte d 17 * 22 ^ . However, the curvatures at small values 
of s (especially as salt concentration decreases) suggests that the tangent correlation function 
qualitatively changes for PEs. The data in Fig. [l]i in a log-log scale clearly shows the postulated 
power law (Eq. Q) at all values of /. The range over which the power law decay is observed is 
from s/a = 1 up to a few A^'s, which is more robust than we conservatively anticipated. Note 
that the fit using Eq. (6a) in ref.^ with a double exponential form, although good for large s 
(Fig. [T]3, solid line), does not reproduce the true behavior in this regime (Fig. [T]i, solid line). 
Thus, the simulations validate Eq. Q. We fit the initial power law decay to extract C and 7 as 
functions of / (or A^)- As shown in Fig. [ij 7 ~ 0.64 (a/A^) ' 36 , whereas C is almost a constant 
for the range of / considered here. 

In order to ascertain the range over which the x ~ In/ behavior is observed using our 
model PE, we performed simulations of a chain under tension using / = 20mM and iV = 
100, 200, . . . , 3200. As shown in Fig. [3^i, the FECs are almost independent of chain length at 
/ ~ 5pN but are clearly iV-dependent for / < 5pN. At even smaller forces, i.e. / < 0.5pN, 
the discrepancies between the N = 1600 and N = 3200 curves are still discernible, which imply 
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strong finite size effects. Nevertheless, the prominent feature of the curves is the x ~ In/ 
behavior which spans almost 3 decades in force (0.1pN< / < 50pN) for all values of N, which 
becomes increasingly transparent as N increases. Fig. |3Jd shows the FECs of a chain with 
N = 3200 at several ionic strengths from ImM to 3000mM. For reference we also show the 
analytic fitP' 



x^l-&/V^ + « 2 (10) 

with a = 2 and c = 4 in the x > 1/2. It can be seen that the x ~ In / regime is visible in 
most of the curves although it is narrower with increasing ionic strength. Moreover, all the 



curves to converge at / ~ 50pN, or x ~ 0.89, which can be analytically fit using Eq. (10). This 
implies that all the chains start to enter the FJC regime and electrostatic interactions among 
the monomers and the ionic conditions play an insignificant role in the high stretching limit. 
Thus, the observed unusual FEC (x ~ In /) is pronounced over a few decades in force /, and is 
surprisingly found to be a consequence of the previously unnoticed power law decay of (cos#(s)) 
(Eq. 4) at short distances. 

Analysis of Experimental Data 

Next we tested the validity of Eq.Q by fitting it to the simulated FECs with only C* and 
c/7 as two adjustable parameters. Because of the approximation made in deriving the equation, 
we restrict the fitting region to be 0.5 < x < 0.8. Indeed, as shown by the solid lines in 
Fig. [3J3 (/ = 20mM) and Fig. [3ji (/ = 50mM), Eq.([8]) fits the data in the region with reasonable 
accuracy. At / = 20mM, the fitted value for c/7 is ~ 4.01 and for C* is ~ 1.78. If we reconcile 
the fitted 7 with the value obtained through independent simulation for the tangent correlation 
at the same I in Fig. [2j we get c = 1.7 and C ~ 0.80 (C is larger than the direct-fit value of 
0.54). Similar values were also obtained at / = 50mM. 

Now we use the theory and simulation results to analyze the stretching data of ssDNA by 
fitting the FECs. Fig. 4 shows the fit of experimental FECs of ssDNA 12 at a few representative 
ionic strengths / using simulated FECs of our PE model, with the contour length L c as the only 
adjustable parameter. The simulation results quantitatively reproduce the experimental FECs. 
The fits at other / from ImM to 3000mM are also excellent. It should be emphasized that the 
sole adjustable parameter L c , which varies from one FEC to another, and is very close to the 
values obtained from the rule from simulated curves at various /, x(f = 50pN) w 0.89 (see 
Fig. 0b). 
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Discussions: 

It should be noted that the Hamiltonian ^ only takes into account the intrinsic bending 
rigidity of the polymer and the screened electrostatic interactions as in the Debye-Huckel the- 
ory The ssDNA molecules used in experiments^ were specially treated to prevent base-paring 
interactions, so it is reasonable not to include them in our model. The high quality of the fits 
using our theory strongly suggests that long-ranged interactions predominantly due to PE ef- 
fects give rise to the x ~ In / behavior . It further validates our arguments that the logarithmic 
dependence of extension on the force is intimately related to the power law decay in the tangent 
correlation at small arc- length separations in the absence of force. Thus, we can readily use the 
values of C and 7 in Fig. [2j that have been extracted from direct simulations of the tangent 
correlation, as predictions characterizing the decay of tangent correlation of ssDNA. 

We can also the use the approximate stretching equation (|7| with c = 1.7 and L c obtained 
from the assumption x(f = 50pN) 0.89, and adjust C and 7 to fit to the experimental FECs. 
Using this procedure we also obtained good agreement with the values obtained from a direct fit 
of Eq. ^ to the simulation data of tangent correlation. The fitted values of C to the experiment 
data are somewhat higher than those obtained using the direct fit, but it also stays almost a 
constant for all the range of I considered. Thus, the analytical theory with only C and 7 as 
parameters suffices to fit the experimental data over a wide range of salt concentration. 

Logarithmic dependence of ssDNA upon stretching has been reported in earlier single 
molecule experiments^. However, in fitting the experimental FECs, they used a different model, 
namely the extensible FJC with electrostatic interactions that are similar to ours in Eq. ([9]). 
As noted here, they also showed that the entropic worm-like chain modeP^ failed to fit the 
experiment data and tended to give unphysical persistence length. However, as pointed out 
elsewhere^, both the discreteness (which is not present in the WLC model, but is in our model 
and the one by Dessinges et alP^) and the bending stiffness play important roles in the elastic 
properties of biomolecules in usual ranges of stretching force. Moreover, we have seen that at 
low ionic strengths, the PE effects in ssDNA are very strong in the range of force less than about 
lOpN (/-dependent). In this range of forces the nature of FECs is mainly determined by the 
electrostatic interactions. At the force higher than lOpN, the semiflexibility of the chain should 
play an important role and the WLC modeP^can be used to quantitatively describe the FECs. 

Finally, it is interesting to wonder why the x ~ In/ dependence has not been observed in 
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dsDNA, which is also an excellent model polyelectrolyte. In dsDNA, the intrinsic persistence 
length (~50nm in monovalent salts) is usually much larger than the Debye length, Xp. Hence, 
the range over which the initial power law decay of the tangent correlation can be observed is 
far too small. For all practical purposes, (cos#(s)) exhibits the standard exponential decay for 
all s. As a result, dsDNA behaves as a WLC with polyelectrolyte effects playing only a minor 
role especially when the salt concentration is higbP^HI Thus, the x ~ In/ behavior cannot be 
observed in practice in monovalent solutions even though in theory such a regime must exist. 
In a recent paper, ssDNA in divalent cations such as Mg 2+ or Ca 2+ shows even a stronger 
dependence of extension x on f than x ~ In / 27 . In order to address these issues, we need to 
develop a better description of ion-PE interactions. 
Conclusions: 

In summary, we have shown that the unusual stretching behavior of ssDNA is due to poly- 
electrolyte effects whose origins can be traced to a previously unnoticed power law decay of the 
tangent correlation function in PEs in the absence of force on length scales comparable to the 
Debye length. Using theoretical arguments and simulations of a minimal model for a PE we 
demonstrated that at forces (10-50 pN) the FEC is given by x ~ In /. Our theory quantitatively 
explains the measured FECM^l (x ~ In/). The simple analytic expression for FECs can be used 
to predict the response of force in other single stranded RNA and DNA as well as synthetic 
polyelectrolytes. 

Acknowledgements: We are grateful to the National Institutes of Health (GM 089685)for 
supporting this research. We thank Omar Saleh for discussions and for providing us with his 
experimental data and a preprint prior to publication. 

Note Added While this paper was under review we received a preprint by M. J. Stevens, 
D. B. Mcintosh, and O. Saleh reporting simulations that illustrate the relevance of crumpled 
structures at short length scales in polyelectrolytes in leading to the observed scaling x ~ In /. 
The physics underlying the origin of x ~ In / in both the treatments are similar. 
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Figure captions: 



Figure 1: (a) A cartoon of the polymer chain showing tangent vectors at s« and Sj along 
the chain. The angle 8(s — \sj — between the two vectors is shown. The average value over 
all possible conformation of the chain (cos#(s)) is the tangent correlation along the chain. The 
length scale over which (cos#(s)) decays at / = is measure of the chain persistence length, 
(b) At / = ImM, all the curves with various N values (from bottom to top N = 100, 200, 800, 
1600) show that there are strong finite size effects. Only when N exceeds 1600 do the curves 
converge, (c) (cos8(s)) in log-linear scale for iV = 1600 and at various /. The solid line is a 
double exponential fit using (cos9) = (1 — (3) exp (— s/Ai) + (3 exp (— s/A 2 ) 24 , with (3 ~ 0.358 
and Ai ~ 1.02a and A 2 ~ 179.11a. (d) To identify the power law behavior the data in (c) is 
shown as a log-log plot. The curves at various ionic strength clearly show that (cos#(s)) decays 
as a power law on scale s < Xd- The arrows correspond to Xd- The solid line is the double 
exponential fit. 

Figure 2: The values of C and 7 extracted from the simulation data of (cos#(s)) for PE 
chains with N = 1600 using Eq. Q. 

Figure 3: (a) Force-extension curves (FECs) for various iV values at / = 20mM. The N 
values from top to bottom are 100, 200, 400, 1600 and 3200 (solid line). The inset shows the 
FEC for N = 3200. (b) FECs for a PE chain with N = 3200 at various values of /. The values 
of / from bottom to top are 1, 10, 20, 50, 100, 200 and 3000mM. The dashed line is plotted using 



x ~ 1 — ^tl yX\ + a 2 Eq. (10). (c) and (d): Fits using Eq. ([8]) to simulation data for a chain 
with N = 3200 at / = 20 and 50mM, respectively. See text for the values of the parameters. 

Figure 4: Single-parameter (using the contour length L c as the fit parameter) fits of sim- 
ulated FECs to experimental data of ssDNA in Na + solutions^. The high quality of the fits 
for / = ImM (a), lOmM (b), lOOmM (c), 500mM and lOOOmM (d) indicates that electrostatic 
interactions play a dominant role in the stretching of ssDNA. 
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